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Abstract 

For an undirected n- vertex planar graph G with non-negative edge- weights, we consider the 
following type of query: given two vertices s and t in G, what is the weight of a min st-cut in G? 
We show how to answer such queries in constant time with 0(nlog^ n) preprocessing time and 
0(n log n) space. We use a Gomory-Hu tree to represent all the pairwise min cuts implicitly. 
Previously, no subquadratic time algorithm was known for this problem. Since all-pairs min cut 
and the minimum cycle basis are dual problems in planar graphs, we also obtain an implicit 
representation of a minimum cycle basis in 0(n log^ n) time and O(nlogn) space and an explicit 
representation with additional 0{C) time and space where C is the size of the basis. 

These results require that shortest paths be unique. Wc deterministically remove this as- 
sumption with an additional log^ n factor in the running time. 
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1 Introduction 



A minimum cycle basis is a minimum-cost representation of all the cycles of a graph whereas the 
all-pairs min cut problem asks to find all the minimum cuts in a graph. In planar graphs the 
problems are intimately related (in fact, equivalent [7]) via planar duality. In this paper, we give 
the first sub-quadratic algorithm for these problems, running in 0(n log^ n) time. In the following, 
we consider undirected graphs with non-negative edge weights. 

All-pairs min cut The all-pairs min cut problem is to find the minimum sf-cut for every pair s, t 
of vertices in a graph G. Gomory and Hu [6j showed that these minimum cuts can be represented 
by an edge- weighted tree such that: 

• the nodes of the tree correspond one-to-one with the vertices of G, 

• for any distinct vertices s and t, the minimum edge-weight on the unique, simple s-to-t path 
in the tree has weight equal to the min st-cut weight in G, and 

• removing the corresponding minimum-weight edge from the tree creates a partition of the 
nodes into two sets that is a partition of the vertices in G corresponding to a min st-cut. 

We call such a tree a Gomory-Hu tree or GH tree; GH trees are also called cut-equivalent and cut 
trees in the literature. Gomory and Hu also showed how to find such a tree with n — 1 calls to a 
minimum cut algorithm. To date, this is the best known method for general graphs and results 
in an 0{n? logn)-time algorithm for planar graphs using the best-known algorithm for min st-cuts 
in planar graphs [3 |9l Hj. There is an algorithm for unweighted, general graphs that beats the 
n — 1 times minimum cut time bound [3j; the corresponding time for planar graphs, however, is 
0(n^poly logn) time. 

Minimum cycle basis A cycle basis of a graph is a set of independent cycles. Viewing a cycle 
as an incidence vector in {0, 1}^, a set of cycles is independent if their vectors are independent 
over GF{2). The weight of a set of cycles is the sum of the weights of the cycles. The minimum- 
cycle basis (MCB) problem is to find a cycle basis of minimum weight. This problem dates to the 
electrical circuit theory of Kirchhoff [13] in 1847 and has been used in the analysis of algorithms 
by Knuth [T5]. For a complete survey, see [10 j . The best known algorithm in general graphs takes 
0{m^) time where oj is the exponent for matrix multiplication |2j. 

The best MCB algorithms for planar graphs use basic facts of planar embeddings: An embedded 
planar graph is a mapping of the vertices to distinct points and edges to non-crossing curves in the 
plane. A face of the embedded planar graph is a maximal open connected set of points that are 
not in the image of any embedded edge or vertex. Exactly one face is unbounded and it is called 
the infinite face. We identify a face with the embedded vertices and edges on its boundary. For 
a simple cycle C in a planar embedded graph G, let int{C) denote the open bounded subset of 
the plane defined by G. Likewise define ext{G) for corresponding unbounded subset. We refer to 
the closure of these sets as int{C) and ext{G), respectively. We say that a pair of faces of G are 
separated by C in G and that G separates this pair if one face is contained in int{G) and the other 
face is contained in ext{C). A set of simple cycles of G is called nested if, for any two distinct 
cycles G and G' in that set, either int{G) C int{G'), int{G') C int{G), or int{G) C ext{G'). 

Hartvigsen and Mardon ^ prove that if G is planar, then there is a minimum cycle basis 
whose cycles are simple and nested in the drawing in the embedding. As such, one can represent a 
minimum cycle basis of a planar embedded graph as an edge-weighted tree, called the MCB tree, 
such that: 



1 



• the nodes of the tree correspond one-to-one with the faces of the planar embedded graph, and 

• each edge in the tree corresponds to a cycle in the basis, namely the cycle that separates the 
faces in the components resulting from removing said edge from the tree. 

Hartvigsen and Mardon also gave an 0(n^ log n)-time algorithm for the problem that was later 
improved to O(n^) by Amaldi et. al. [2]. 

1.1 Planar duality 

In planar graphs, the MCB and GH problems are related via planar duality. 

Corresponding to every connected planar embedded graph G (the primal) there is another 
connected planar embedded graph (the dual) denoted G*. The faces of G are the vertices of G* 
and vice versa. Two vertices in G* are connected if the corresponding faces in G share a border. 
Note that the edges of G correspond one-to-one with those of G* . Cycles and cuts are equivalent 
through duality: 

In a connected planar graph, a set of edges forms a cycle in the primal iff it forms a cut 
in the dual. [H] 

Equivalence between MCB and GH trees Just as cuts and cycles are intimately related via 
planar duality, so are the all-pairs min cut and minimum cycle basis problems. In fact, Hartvigsen 
and Mardon showed that they are equivalent in the following sense: 

Theorem 1. For a planar embedded graph G, a tree T represents a minimum cycle basis of G if 
and only if T is a Gomory-Hu tree for G* (after mapping the node-face relationship to a node-vertex 
relationship via planar duality). (Corollary 2.2 'Jj) 

Herein, we focus on the frame of reference of the minimum cycle basis. Our algorithm works by 
finding a minimum (weight) cycle that separates two as- yet unseparated faces / and g. By duality, 
this cycle is a min fg-cut in G* and for conciseness we name separating cycles in G min-cuts for 
the respective faces. We recurse on unseparated faces, gradually building the tree T that represents 
the minimum cycle basis. By Theorem [T| T is the GH tree for the dual graph. This alone will not 
achieve a sub-quadratic running time. In order to beat quadratic time, we guide the recursion with 
planar separators and use precomputed distances to efficiently find the minimum separating cycles. 

1.2 Planar separators 

A decomposition of a graph G is a set of subgraphs Pi, . . . ,Pk such that the union of vertex sets 
of these subgraphs is the vertex set of G and such that every edge of G is contained in a unique 
subgraph. We call Pi, . . . ,Pk the pieces of the decomposition. The boundary vertices of a piece Pi 
is the set of vertices u in that piece such that there exists an edge {u, v) in G with v ^ Pi. We use a 
recursive subdivision. When piece P is decomposed into child subpieces, the boundary nodes of a 
child are the boundary vertices inherited from P as well as the boundary vertices introduced by the 
decomposition of P. By recursively applying Miller's Cycle Separator Theorem [16j as described 
in [5] until each leaf piece has a single edge, one obtains the following: 

Definition 1 (Balanced recursive subdivision). A decomposition of G such that at each level a 
piece with n nodes and r boundary nodes is divided into two subpieces each containing no more 
than |n -|- c^/n nodes and no more than |r + c^/n boundary nodes. 
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Fixing an embedding of G, a piece inherits an embedding from G's embedding. A hole in a 
piece is a bounded face containing boundary vertices. While it is not possible to guarantee that 
holes are not introduced by the balanced recursive subdivision, for the simplicity of presentation, 
we assume that such a holeless decomposition exists. Hence, all boundary vertices of a piece are 
on its external face. We will also assume that pieces are connected. We show how to handle the 
general case in Section [7j 

We define the O(logn) levels of the recursive subdivision in the natural way: level consists of 
one piece (G) and level z-pieces are obtained by applying the Cycle Separator Theorem to each level 
i — 1-piece. We represent the recursive subdivision as a binary tree, called the subdivision tree (ofG), 
with level i-pieces corresponding to vertices at level i in tree. Parent/child and ancestor /descendant 
relationships between pieces correspond to their relationships in the subdivision tree. 

Dense distance graphs For a piece P, the internal dense distance graph of P or miDDG(P) 
is the complete graph on the set of boundary vertices of P, where the weight of each edge (u, v) 
is equal to the shortest path distance between u and v in P. The union of internal dense distance 
graphs of all pieces in the recursive subdivision of G is the internal dense distance graph (of G), or 
simply iniDDG. Fakcharoenphol and Rao showed how to compute intDDG in 0(n log^ n) time [5j; 
Klein improved this to O(Tilog^n) |14) . 

The external dense distance graph of P or extDI}G{P) is the complete graph on the set of 
boundary vertices of P where the weight of an edge (n, v) is the shortest path distance between u 
and V m. G\ E{P). The external dense distance graph of G or exiDDG is the union of all external 
dense distance graphs of the pieces in the recursive subdivision of G. 

Theorem 2. The external dense distance graph of G can he computed in O(nlog^n) time. 

Proof. Fakcharoenphol and Rao compute intDDG bottom-up by applying a variant of Dijkstra to 
obtain iniDDG(P) for a piece P from the internal dense distance graphs of its two children. A 
similar algorithm can compute extDDG(P) from the external dense distance graph of the parent of 
P and the internal dense distance graph of the sibling of P. Hence, we can obtain exiDDG with a 
top-down algorithm after having found in^DDG. The running time is the same as that for finding 
mtDDG. □ 

1.3 Overview of the algorithm 

We gradually build up the tree representing the minimum cycle basis. Initially the tree is a star 
centered at a root r and each leaf corresponding to a face in the graph (including the infinite face). 
We update the tree to reflect the cycles that we add iteratively to the basis. When the first cycle 
C is found, we create a new node xc for the tree, make xc a child of the root and make all the 
faces that C encloses children of xc- Each non-face node in the tree corresponding to a cycle C 
defines a region R: the subgraph of G contained in the closed subset of the plane defined by the 
interior of C and the exterior of the children (if any) of G. We say that R is hounded by C, that G 
is a hounding cycle of R, and that R contains the child regions and/or child faces defined by the 
tree. The tree of regions is called the region tree. Observe that a pair of faces not yet separated by 
a basis cycle belong to the same region. 

The root r will remain a special region that represents the entire plane. We only add cycles to 
the basis that nest with the cycles found so far. Whenever the basis is updated, the region tree is 
updated accordingly. This is illustrated in Figure [Tj We show how to efficiently update the region 
tree in Section [H 
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In the final tree, all faces have been separated: each face is the only face-child of a region. We 
call such a region tree a complete region tree. Mapping each face to its unique parent, creating a 
tree with one node for each face in the graph, will create the MCB-tree. 
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Figure 1: A graph with faces a through g; four nesting cycles A through D (left). A region tree for 
cycles A,B and C (center). A region tree for cycles A through D (right). 

Our algorithm is guided by the recursive subdivision of G. Starting at the deepest level of the 
recursive subdivision, we separate all pairs of faces of G that have an edge in a common piece of 
the subdivision. Each (nesting) separating cycle is added to the region tree. 

Suppose we are at level i and that our algorithm has been applied to all levels deeper than i: 
every pair of faces that have an edge in a common subpiece have been separated. Consider some 
level i-piece P and suppose two faces /i and /2 of G both sharing edges with P have not yet been 
separated. Let Pi and P2 be the two child pieces of P in the recursive subdivision. Since all pairs 
of faces sharing edges with Pi and all pairs of elementary faces sharing edges with P2 have already 
been separated, w.l.o.g. /i shares edges with Pi and /2 shares edges with P2. 

Since /i and /2 have not been separated, they must belong to a common region R in the region 
tree. There cannot be three unseparated faces of R in P. For otherwise assume w.l.o.g. that there 
is a third unseparated face / that shares edges with Pi . Since all pairs of faces sharing edges with 
Pi have already been separated and since both fi and / belong to P, / = /i a contradiction. 

Let the region subpieces of a piece P be the subgraphs defined by the non-empty intersections 
between P and regions defined by the region tree. We say that a region R is associated with a 
region subpiece Pr and that Pr is associated with R if Pr = P n P is not empty. The boundary 
nodes OPr of Pr are the boundary nodes that are inherited from P. The above shows that in 
each region subpiece Pr of P, at most one pair of faces needs to be separated. These constructs 
are illustrated in Figure [2] In Section [3j we show how to separate such a pair of faces in time 
0(|Pr| log^ \Pr\) assuming when two subpieces remain connected. When this assumption does not 
hold, we need 0(|Pr| log'^ \ Pr\) time to handle the general case as shown in Section [rj This time 
amounts to 0(|P| log^ |P|) time over all region subpieces of P. This will imply that our algorithm 
spends O(nlog'^n) time at level i in the recursive subdivision. Summing over all levels, this gives 
a total running time of 0(n log^ n). 



1.4 Results 

The bulk of the paper will be devoted to presenting the details of the region tree algorithm outlined 
in the previous section. The main technical part of the algorithm is presented in Section [3| where 
we show how to find the minimum separating cycle for each of those pairs of unseparated faces in 
such a way that the total time over all pieces is bounded by 0{n). However, this assumes that we 
can efficiently find all region subpieces that contain unseparated faces. We describe how to do so 
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Figure 2: The dotted edges are the edges belonging to the boundaries of a piece P and P's children 
Pi and P2. f is a face of Pi and g is a face of P2. The solid black edges are the bounding cycle of 
a region R. The three shaded regions are three child regions of R. The thick grey edges are the 
edges of region subpiece Pr. (The remaining edges of the graph are the thin, grey edges.) The 
intersection of a minimum separating cycle for / and g with P uses only edges of Pr. 

in Section |4j In Section [5} we show how to update the region tree given a new nesting separating 
cycle. In Section [7| we show how to handle general case when holes are present, and when a piece 
is decomposed into more then two connected subpieces. Together these four sections will prove the 
main result of our paper: 

Theorem 3. A complete region tree of a planar undirected n-vertex graph with non-negative edge 
weights can he found in O(nlog^n) time and 0(n log n) space. 

As argued above the complete region tree can be used to obtain the MCB-tree. Therefore, we 

get: 

Theorem 4. The minimum cycle basis of an undirected and unweighted planar graph can be com- 
puted in O(nlog^n) time and 0(n log n) space. 

To compute the Gomory-Hu tree of a connected planar graph G, find the MCB-tree T of G*. 
By Theorem [T| mapping faces of G* represented in T to their corresponding vertices in G produces 
the GH tree of G, giving: 

Theorem 5. A Gomory-Hu tree of an n-vertex connected undirected planar graph with non-negative 
edge weights can be computed in O(nlog^n) time and 0(n log n) space. 

In order to find a minimum st-cut using the GH tree, we need to find the minimum weight edge 
on the s-to-t path in the tree. With an additional 0(n log n) preprocessing time, one can answer 
such queries in 0(1) time using a tree-product data structure [12], giving: 

Theorem 6. With O(nlog^n) time and 0(n log n) space for preprocessing, the weight of a min 
st-cut between for any two given vertices s and t of an n-vertex planar, undirected graph with 
non-negative edge weights can be reported in constant time. 

In Section |6j we will show how to explicitly find the cycles given the complete region tree, giving 
the following results: 
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Theorem 7. Without an increase in preprocessing time or space, the min st-cut oracle of Theorem^ 
can he extended to report cuts in time proportional to their size. 

Theorem 8. The minimum cycle basis of an undirected planar graph with non-negative edge weights 
can be computed in 0{n log^ n + C) time and 0{n log n + C) space, where C is the total size of the 
cycles in the basis. 

In the last section (Section [s]) we derandomize the algorithm by paying an additional log^ n 
factor in the running time. 



2 Preliminaries 

To simplify the presentation of the algorithm and the analysis, we make a few structural assumptions 
on the input graph. These assumptions are not truly restrictive. 



2.1 Simplifying structural assumptions 

Triangulation We assume that a the input we are given a triangulation, i.e., a planar graph in 
which all faces are triangles. We can always achieve this by triangulating the graph with infinite- 
weight edges. For each face / of the original graph, identify a face /' of the triangulated graph that 
is enclosed by / in the inherited embedding. The min f'g'-cut will not use any infinite-weight edge: 
the cut contains the same set of edges as in the min fg-cut in the original graph. Likewise, the 
set of cycles in a minimum cycle basis in the original graph are mapped to the set of finite-weight 
cycles in a minimum-cycle basis of the triangulated graph. As a result of the triangulation we get 
that each face in the graph is simple: each vertex appears only once on the (boundary of the) face. 



Degree three We assume that each vertex has degree three. This can be achieved by triangulat- 
ing the faces of the dual graph, giving a degree three primal graph. We triangulate the dual graph 
with e weight edges, where e is much smaller than the smallest weight 6e of an edge e.g., e = ^ 
times smaller. We use e weight rather than weight to aid in the next structural assumption. As 
in the above triangulation, we can map between minimum cycle bases and min cuts in the original 
graph and the degree-three graph. Triangulating the dual will increase the size of the faces of the 
primal, but the faces will remain simple: assuming that faces are simple and that vertices have 
degree three are compatible assumptions. Each vertex v in the original graph is mapped to a tree 
Ty of e-weight edges in the degree-three graph. 



Unique shortest paths We assume that in the graphs we consider (G as well as the cut-open 
graphs in Section |3]), there is a unique shortest path between any pair of vertices. To make 
this assumption a reality one can add a small, random perturbation to the weight of each edge 
and make the probability of having non-unique shortest paths arbitrarily small. For example, if 
we take a random perturbation from the set [1, . . . ,n^]^, then by the Isolating Lemma due to 
Mulmuley, Vazirani and Vazirani |17j the probability that all shortest paths are unique is at least 
1 — Kx (T\ > 1 — -. Moreover, since this random fraction is much smaller than the e used above, we 
guarantee that any shortest path that enters and leaves a tree of e-weight edges does so exactly 
once. Theorem 11 in Section [8] gives a more robust, deterministic way to avoid this assumption 
that requires a log^ n-increase in the running time of our algorithm. 
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Every region is bounded by a cycle The root of the region tree described in Section 1.3 was 
a special region corresponding to the entire graph. We can treat this region as any other region by 
bounding the graph by a zero-weight cycle, and so treat every region as having a bounding cycle. 



2.2 Isometric cycles 

In a planar embedded graph, a simple cycle C is said to cross another simple cycle C if the cycles 
are not nested. 

A cycle C in a graph is said to be isometric if for any two vertices u,v E C, there is a shortest 
path in the graph between u and v which is contained in C. A set of cycles is said to be isometric 
if all cycles in the set are isometric. The following two results will prove useful. 

Lemma 1. Assume that shortest paths in a graph G are unique. The intersection between an 
isometric cycle and a shortest path in G is a (possibly empty) shortest path. The intersection between 
two distinct isometric cycles C and C' in G is a (possibly empty) shortest path; in particular, if G 
is a planar embedded graph, C and C' do not cross. 

Proof. Let C be an isometric cycle and let P be a shortest path intersecting C. Let u resp. v be 
the first resp. last vertex of intersection between C and P for some orientation of P. Since C is 
isometric, there is a shortest path P' in G between u and v. Since shortest paths are unique, P' is 
the subpath of P between u and v. Hence, the intersection between C and P is shortest path P' . 
This shows the first part of the lemma. 

Let C" 7^ C be an isometric cycle. Then C H C is a set of vertex-disjoint paths. We claim that 
there can be at most one such path. To see this, let u and v be distinct vertices on paths P„ and 
Pv Let Pi and P2 be the two edge-disjoint paths in C between u and v. Define P[ and P2 similarly 
for C' . Since C and G' are isometric, we may assume that, say. Pi and P[ are shortest paths. By 
uniqueness of shortest paths. Pi = P{ so Pu = Pv, implying that C Pi C" is a path. If we pick u and 
V as the first and last vertex on this path, it follows that C H C" is a shortest path. □ 

Lemma 2. Any minimum cycle basis of a graph is isometric. (Proposition 4-4 ^) 

It follows from these two lemmas that the minimum cycle basis that our algorithm constructs 
is isometric and nested. 

2.3 Representing Regions 

We will represent the region tree using the top tree data structure jl]. This will allow us to find 
the lowest common ancestor lca{x, y) of two vertices x and y and determine whether one vertex is 
a descendant of another in logarithmic time and also to update the region tree efficiently as we add 
cycles to the basis (the latter is described in Section [sj . We will also make use of the operation 
jump{x,y,d), which for two vertices x and y in a top tree finds in logarithmic time the vertex of 
distance d (in terms of edges) from x on the simple path between x and y. 

If the top tree represents a weighted tree, it can report the weight of a simple path in logarithmic 
time, given the endpoints of this path. 

Additionally, for each region in the region tree we store its compact representation. A compact 
representation of a region is G[P], where vertices with degree 2 are removed by merging the ad- 
jacent edges creating super edges. For each super edge we store the first and the last edge on the 
corresponding path. We show how to maintain compact representations of regions in Section [5] 
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3 Separating a pair of faces 



In this section we show how to find the minimum /(^-separating cycle for the unique pair of un- 
separated faces /, g in a region subpiece Pr. We assume that the region subpiece is given to us by 
the work in Section [4] We emulate the algorithm due to Reif |18j to find the minimum separating 
cycle. 

We will use an operation of cutting open a planar embedded graph G along a path X: duplicate 
every edge of X and every internal vertex of X and create a new, simple face whose boundary is 
composed of edges of X and their duplicates. The resulting graph is denoted Gx- 

Paths P and Q cross if there is a quadruple of faces adjacent to P and Q that cover the set 
product {left of P, right of P] x {left of Q, right of Q}. 



3.1 Reif's minimum separating cycle algorithm 

Let X be the shortest path between any vertex on the boundary of / and any vertex on the boundary 
of g. Reif uses the fact that, since X is a shortest path, there is a minimum (j(/-separating cycle, 
C, that crosses X only once: 

Theorem 9. Let X he the shortest g-to-f path. For each vertex x € X , let Gx be the minimum 
weight cycle that crosses X exactly once and does so at x. Then a minimum fg-separating cycle 
is a cycle Gx of minimum weight. Further, Gx is the shortest path between duplicates of x in Gx ■ 
(Proposition 3 /7^ . originally due to Itai and Shiloach /71]/. ) 

Reif computes, for every vertex x G X, the minimum separating cycle Gx- Finding Gx amounts 
to finding a shortest x-to-x' path in Gx, where x' is the duplicate of x. The running time of the 
algorithm is bounded by divide and conquer: start with x, the midpoint of X in terms of the 
number of vertices, and recurse on the subgraphs obtained by cutting along Gx- The algorithm 
takes O(nlogn) time using the linear-time algorithm for shortest paths in planar graphs. 

One of the main ingredients of Fakcharoenphol and Rao (Section 3.2.2 and 3.2.2 [5]) is the im- 
plementation of Dijkstra algorithm that executed on a piece P runs in timeO((|5P| + |P|)log2|P|). 
This allows us to emulate Reif's algorithm in time 0((|5Pr| + \Pr\) log^ I-PrD- In order to attain 
this running time, we must deal with the following peculiarities: 

• X is not contained entirely in Pr. For a vertex x £ X that is outside Pr, we will compute 
Gx by composing distances in exiDDG and intDDG between restricted pairs of boundary 
vertices of Pr. We call such cycles external cycles. We show how to find these cycles in 
Section 13. 4[ 



For vertices x G X that are inside Pr, Gx may not be contained by Pr. We find Gx by first 
modifying extDDG to disallow paths from crossing X. We call such cycles internal cycles. 



We show how to find these cycles in Section 3.3 



extDDG corresponds to distance in G, not Gx- We compute modified dense distance graphs 
to account for this in Section [3.21 



We can compute X with using Dijkstra in 0((|9Pr| + |Pr|) log^ \ Pr\) time. In the next two sections 
we show how to find all the internal cycles and external cycles in time 0{{\dPR\ + \Pr\) log^ \Pr\) 
time. The minimum length cycle over all internal and external cycles is the minimum /(7-separating 
cycle in G. 
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Figure 3: An external and internal cycle separating faces / and g in a. region subpiece, whose 
boundary vertices are given by the bold vertices. 

3.2 Modifying the external dense distance graph 

We use extDDG to compute extDDGx, the dense distance graph that corresponds to distances 
between boundary vertices of Pr when the graph is cut open along X. However, we do not 
compute extDDGx as it can take too much time. We simply show how to determine its values 
when needed. 

Let B be the set of boundary vertices of Pr. Cutting G open along X duplicates vertices of 
B that are in X, creating B' . extDDGx can be represented as a table of distances between every 
pair of vertices of B': extDDGx (a;, y) = oo if x is a copy of y or if x and y are separated in Gx 
outside Pr and extDDGx ix,y) = extDDG{x,y) otherwise. 

We describe how to determine if x and y are separated in Gx outside Pr. The portions of X 
that appear outside Pr form a parenthesis of (a subset of) the boundary vertices, illustrated in 
Figure |4j By travelling along X we can label the start and endpoints of these parentheses. By 
walking along the boundary of the subpiece we can label a group of boundary vertices that are not 
separated by X by pushing the vertices onto a stack with a label corresponding to the start of a 
parenthesis and popping them off when the end of the parenthesis is reached, labelling the boundary 
vertices with the corresponding parenthesis. Two boundary vertices are not separated if they have 
the same parenthesis label. Hence, whenever we are asked for a distance in extDDG x{x,y) we 
either answer oo if x and y are not in the same parenthesis, or return extDDG{x, y) otherwise. 




Figure 4: Modifying the external dense distance graph. (Left) X is given by the solid line and the 
boundary of the subpiece is given by the dotted line. The parts of X outside the subpiece form a 
parenthesis. (Right) In Gx, the only finite distances from a in extDDGx correspond to the thick 
lines. The shaded area represents the new face created by cutting along X. 
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3.3 Finding internal cycles 



Let D = V{X n Pr). Let D be ordered according to the order of the vertices along X. For 
each vertex x £ D, we compute the shortest x-to-x' paths in Gx where x' is the copy of x in 
Gx- We do this using Dijsktra's algorithm on the cut-open graph induced by the vertices in 
Pr (i.e. Gx[Pr]) and the modified dense distance graph: extDDGx- Each cycle can then be 
found in 0{{\Pr\ + |9PrP) log |Pr|) time. Let Xm be the midpoint vertex of D according to the 
order inherited from X. Gx^ splits Pr and extDDGx into two parts (not necessarily balanced). 
Recursively finding the cycles through the midpoint in each graph parlj^ results in log|D| levels 
for a total of 0{{\Pr\ + I^P/jp) log |Pr| log \D\) = 0{{\Pr\ + I^P/jp) log^ |P^|) time to find all the 
internal cycles. 



3.4 Finding external cycles 

We show below that every external cycle is composed of a single edge ab in the unmodified extDDG 
and a shortest path -Kab between boundary vertices of Pr in G that does not cross X. 

Using the modified Dijkstra algorithm of Fakcharoenphol and Rao (Section 3.2.2 and 3.2.2 [5]), 
all shortest paths TTab can be computed in a total of 0{\dPR\^log^\PR\) time if extBBGx and 
intDDGx are given: for each a, find a shortest path tree in extDT>Gx U intDDGx rooted at a in 
0{\dPR\ log^ \Pr\) time. intDBGx can be found in 0((|9PkP + I-PrI) log^ I^'hI) time by cutting 
open X and using the algorithm of Fakcharoenphol and Rao. Note that we need to compute 
intT>T)Gx from scratch because X has been cut open and because Pr is no longer a subgraph of 



G due to the compact representation presented in Section |4.3 

In order to compute all the external cycles, one enumerates over all pairs a, 6 of vertices that 
are split exactly once by the parenthesis given by X, summing the weight of TTab and the weight of 
edge (o, 6) in extDDG. Since there are 0(|9Pr|) boundary vertices, there are 0(|9PrP) pairs to 
consider. The minimum-weight external cycle then corresponds to the pair with minimum weight. 
By the above, this cycle can be found in ©((I^PrP + |Ph|) log^ I-PrJ) time. 

It remains to prove the required structure of the external cycles. 

Lemma 3. The shortest external cycle is composed of a single edge ab in the unmodified extDDG 
and a shortest path TTab between boundary vertices of Pr in G that does not cross X. 

Proof. Let G be the shortest external cycle that separates faces / and g. By Theorem [9| C is 
a cycle that crosses X exactly once, say at vertex x. Further, C is a shortest path P between 
duplicates of x in the graph Gx- Since G must separate / and g, G must enter Pr. Starting at x 
and walking along G in either direction from X, let a and b be the first boundary vertices that G 
reaches. Let TTab be the a-to-6 subpath of G that does not cross X. Since C is a shortest path in 
Gx, T^ab is the shortest path between boundary vertices as given in the theorem. 

Let TTx be the a-to-6 subpath of G that does cross X. By definition of a and 6, tt^ contains 
no vertices of Pr except a and b. Further, tt^ must be the shortest such path, as otherwise, G 
would not be the shortest /^-separating cycle. Therefore, tTx must correspond to the edge ab in 
extDDG. □ 

In Section [7] we will present a different method for computing external cycles that generalizes 
to disconnected pieces. 



^In order to properly bound the running time, one must avoid repeating long paths in the subproblems: in a 
subproblem resulting from divide and conquer, we remove degree-two vertices by merging the adjacent edges. 
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4 Finding region subpieces 



In Section |1.3[ we defined the region subpieces of a piece as the intersection between a region and 
a piece (Figure [2]). In this section, we show how to identify the region subpieces and the edges that 
are in them. We start by identifying the set of regions TZp whose corresponding region subpieces of 
piece P each contain a pair of unseparated faces (Section 4.1 ). For each region R G TZp we initiahze 
the corresponding region subpiece Pp as an empty graph. For each edge e of P we determine to 
what region subpieces e belongs using lowest common-ancestor and ancestor-descendent queries in 
the region tree (Section 4.2). In Section 4.3 we show how to do all this in 0(r log^ n] 
is the size of P. 



time where r 



4.1 Identifying region subpieces 

Since each edge is on the boundary of two faces, we start by marking all the faces of G that share 
edges with P in 0(r) time. Since a pair of unseparated faces in P are siblings in the region tree, we 
can easily determine the regions that contain unseparated faces. So, in 0(r) time we can identify 
TZp, the set of regions with unseparated faces in P. We will need the following bound on the size 
of T^P in our analysis, illustrated in Figure [5} 




Figure 5: A piece P is given by the boundaries (dotted) by its two child pieces Pi and P2- TZp is 
a nesting set (solid cycles), each containing a pair of unseparated faces (grey). Since these faces 
must be separated by the child pieces, each bounding cycle (except for one outer cycle) in T^p must 
cross the dotted lines, resulting in a bound on \TZp\. 



Lemma 4. \TZp\ = 0{y/r). 

Proof. Let B be the set of boundary nodes in Pi and P2. Let J-i be the set of faces containing 
edges of Pi and not edges of P2 and let F2 be the set of faces containing edges of P2 and not edges 
of Pi. 

Any region R G TZp must contain at least one face of each and J-2- So, if C is the cycle 
bounding P, either int{C) contains B or C crosses B, i.e., crosses the bounding faces of Pi or P2. 

If int{C) contains B then we claim that no other cycle bounding a region in T^p has this property. 
To see this, suppose for the sake of contradiction that there is another such cycle C bounding a 
region R' G TZp. The cycles have to nest, i.e., either int{C) C int{C') or int{C') C int{C). Assume 
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w.l.o.g. the former. Then all faces of Ti are contained in the interior of a face of R' . But this 
contradicts the assumption that R' contains at least one face from Ti. 

We may therefore restrict our attention to regions R^TZp whose bounding cycle C crosses B. 

Consider two vertices u and u in B we and let T^-p" be the set of regions whose bounding cycles 
go through u and v. By Lemma [T] we know that if two isometric cycles share two vertices, then 
they share a path between these vertices as well. Hence, all cycles bounding regions in TV^^ share 
a path vr from u io v. As a result at most two of the regions in T^p*" can contain a face adjacent to 

TT. 

By Lemma [T] no bounding cycles of regions in T^p cross, so each cycle going through non- 
consecutive boundary vertices splits the set of boundary vertices into two parts - one inside and 
the other one outside. By a "chocolate-breaking" argument there cannot be more than — 1 such 
pairs of non-consecutive vertices used by cycles. Moreover, there are no more than \B\ consecutive 
pairs possible. As argued above each of these pairs cannot be used by more than two regions that 
contain a vertex inside P, so there are no more than 4|i?| = 0{\/r) regions in IZp. 

□ 



4.2 Identifying edges of region subpieces 

Region subpieces are composed of two types of edges: internal edges and boundary edges. Let R be 
a region and let C be the bounding cycle of R. An edge e is an internal edge of a region subpiece 
R if the faces on either side of e are enclosed by C. An edge e is a boundary edge of -R if e is an 
edge of C. Every edge is an internal edge for exactly one region subpiece: Lemma [5] shows how 
we can identify this region. We can also determine if an edge is a boundary edge for some region 
(Lemma [6|. However, an edge can be a boundary edge for several region subpieces. We show how 



to deal with this potential problem efficiently in Section 4.3 



Lemma 5. Let e he an edge of G and let fi and f2 be the faces adjacent to e. Then e is an internal 
edge of a region R iff R is the lowest common ancestor of f\ and fi in the region tree. 

Proof. There must exist some region R satisfying the lemma. Let C be its bounding cycle. Then 
both /i and /2 are contained in int{C) and it follows that R must be a common ancestor of /i and 
/2. If R' is another common ancestor then either R is an ancestor of R' or R' is an ancestor of 
R. If R is an ancestor of R' then R' is contained in a face of i? so e cannot belong to R. But this 
contradicts the choice of R. Hence, R = lca{fi, /2). □ 

Iterating over each edge e of P, we can find the region R for which e is an internal edge. If 
R £ TZp, we add e to the corresponding region subpiece Pr. The total time to add internal edges 
to their corresponding region subpieces is 0(r log n). 

Lemma 6. Let e be an edge of G and let fi and f2 be the faces adjacent to e. Let R' be the lowest 
common ancestor of fi and /2 in the region tree. Then e is a boundary edge of a region R iff R is 
a descendant of R' and one o//i,/2 is a descendant of R. 

Proof. Assume first that e G C, where C is the cycle bounding R. Then exactly one of the faces /i 
and /2 is in int{G) and the other is in extiC). Assume w.l.o.g. that fi G int{G) and /2 G ext[C). 
Then /i is a descendant of R and since /2 is not, R must be a descendant of i?'. 

Now assume that i? is a descendant of i?' and that, say, /i is a descendant of R. Then /2 is not 
a descendant of R since otherwise, R' could not be an ancestor of R. This implies that e G C. □ 

Let i? be a region in T^p and let G be the bounding cycle of R. Let Pi and P2 be the children 
of P and let B be the union of boundary vertices of Pi and P2. We have seen (in the proof of 
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Lemma [4]) that the intersection between C and P are subpaths in P between pairs of vertices of 
B. We shah refer to these as cycle paths (of C). Consider the fohowing algorithm to find starting 
points of cycle paths. 

Cycle path starting points identification algorithm Pick a boundary vertex u € B. For 
every edge e adjacent to u (there are at most three such edges), check to see if e is a boundary 
edge of R. If there is no such edge, then there is no cycle path through u. Otherwise, mark 
e as a staring point of a cycle path for R. Repeat this process for every vertex in B. 

Using Lemma [g] and the top tree, this process takes O(y^logn) time for each region R since a 
constant number of tree queries for every vertex in B. By Lemma |4] doing this for all regions in 
TZp takes 0(r log n) time. 

After identifying staring points for cycle-path we can find all edges belonging to them using 
linear search, i.e., the next edge on the cycle C is found by looking an endpoint of the previous 
edge and checking which of the two remaining edges belongs to C. If the cycles are edge-disjoint 
over all regions R € TZp, then the cycle paths will also be edge-disjoint. In such a case the time to 
find all the region subpieces using linear search is 0{{\^lp\^/r + |P|)logn) = O(rlogn). However, 
the cycles are not necessarily edge disjoint. We overcome this complication in the next section. 

4.3 Efficiently identifying boundaries of region subpieces 

Since cycles will share edges, the total length of cycle paths over all cycles can be as large as 
0(r^/^). However, we can maintain the efficiency of the cycle path identification algorithm by 
using a compact representation of each cycle path. The compact representation consists of edges 
of P and cycle edges that represent paths in P shared by multiple cycle paths. 

View each edge of G as two oppositely directed darts and view the cycle bounding a region as 
a clockwise cycle of darts. The following is a corollary of Lemma [T] 

Corollary 1. If two isometric cycles C and C of G share a dart, then either int{G) C int{G') or 
mi{G') C i^(C). 

Let T be the forest representing the ancestor/descendant relationship between the bounding 
cycles of regions in Tip. By Lemma [4| there are 0{y/r) bounding cycles and since we can make 
descendent queries in the region tree in O(logn) time per query, we can find J- in 0(r log n) time. 
Let d be the maximum depth of a node in (roots have depth 0). For z = 0, . . . , d, let Cj be the 
set of cycles corresponding to nodes at depth i in T. As a result of Lemma [T| 

Corollary 2. For any i G {0, . . . ,d}, every pair of cycles in Ci are pairwise dart-disjoint. 
Bottom-up algorithm 

We find cycle paths for cycles in Cd, then C^-i, and so on. The cycles in Cd are dart disjoint, so 
any edge appears in at most two cycles of Q. We find the corresponding cycle paths using the 
cycle path identification algorithm in near- linear time. While Corollory [2] ensures that the cycles 
in Cd are mutually dart-disjoint, they can share darts with cycles in Cd-i- In order to efficiently 
walk along subpaths of cycle paths Q that we have already discovered, we use a balanced binary 
search tree (BBST) to represent Q. We augment the BBST to store in each node the length of the 
subpath it represents. Now, given two nodes in Q, the weight of the corresponding subpath of Q 
can be determined in logarithmic time. 

To find the cycle paths of a cycle C € Cd-i that bounds a region R, we emulate the cycle path 
identification algorithm: start walking along a cycle path Q of G, starting from a vertex of B, and 
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binary 




Figure 6: Finding a cycle path (highlighted straight line) for a cycle C G Q-i between boundary 
nodes of Pi and P2 (grey dashed lines) is found by alternating linear (solid) and binary (dotted) 
searches. Binary searches corrrespond to cycle paths of region subpieces (shaded) bounded by 
cycles in Cd- 

stop if you reach an edge e = uv that has already been visited (linear search). In this case, e must 
be an edge of a cycle path Q' of a cycle C G Cd- By Lemma [T} the intersection of Q and Q' is a 
single subpath and so we can use the BBST to find the last vertex w common to Q and Q' (binary 
search). We add to Pr an edge uw of weight equal to the weight of the u-to-w subpath of Q to 
compactly represent this subpath. If w £ B, we stop our walk along Q. Otherwise we continue 
walking (and adding edges to the corresponding region subpiece) in a linear fashion, alternating 
between linear and binary searches until a boundary vertex is reached. See Figure [6j 

We have shown how to obtain region subpieces for cycles in Cd and in Cd-i- In order to repeat 
the above idea to find cycle paths for cycles in Q_2) we need to build BBSTs for cycle paths of 
cycles in Let Q be one such cycle path. Q can be decomposed into subpaths QiQ'i • • • QkQ'ki 

where Qi, . . . , Qk are paths obtained with linear searches and Q'^, . . . , Q'^. are paths obtained with 
binary searches (possibly Qi and/or Q'^ are empty). To obtain a binary search tree T for Q, we 
start with T the BBST for Qi. We extract a BBST for Q'^ from the BBST we used to find Q'l and 
merge it into T. We continue merging with BBSTs representing the remaining subpaths. 

Once BBSTs have been obtained for cycle paths arising from Cd-i, we repeat the process for 
cycles in Q_2, . . • ,Co. 

Running time We now show that the bottom-up algorithm runs in 0(r log^n) time over all 
region subpieces. We have already described how to identify boundary vertices that are starting 
points of cycle paths in 0(r log n) time. It only remains to bound the time required for linear and 
binary searches and BBST construction. 

A subpath identified by a linear search consists only of edges that have not yet been discovered. 
Since each step of a linear search takes 0(log n) time, the total time for linear searches is 0(r log n). 

The number of cycle paths corresponding to a cycle C is bounded by the number of boundary 
vertices, 0{y/r). We consider three types of cycles paths. Those where 

1. all edges are shared by a single child of C in J^, 

2. no edges are shared by a child, and 

3. some but not all edges are shared by a single child. 

Cycle paths of the first type are identified in a single binary search which, by Lemma [4j sums up 
to a total of 0{r) binary searches over all cycles C £ T . Cycle paths of the second type do not 
require binary search. For a cycle path Q in the third group, Q can only share one subpath with 
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each child (in F) cycle by Lemma [T| hence, there can be at most two binary searches per child. 
Summing over all such cycles, the total number of binary searches is 0{\/r) by Lemma [4] 

In total there are 0{r) binary searches. Each BBST has 0{r) nodes. In traversing the binary 
search tree, an edge is checked for membership in a given cycle path using Lemma |6] in O(logn) 
time. Each binary search therefore takes O(logrlogn) = O(log^n) time. The total time spent 
performing binary searches is O(rlog^n). 

It remains to bound the time needed to construct all BBSTs. We merge BBSTs Ti and T2 in 
0(min{|ri|, IT2I} logdTil + IT2I})) = 0(min{|Ti|, IT2I} logn) time by inserting elements from the 
smaller tree into the larger. 

When forming a BBST for a cycle path of a cycle C, it may be necessary to delete parts of 
cycle paths of children of C. By Lemma [T| these parts intersect mt{C) and will not be needed for 
the remainder of the algorithm. The total number of deletions is 0(r) and they take 0(r log r) 
time to execute. So, ignoring deletions, notice that paths represented by BBSTs are pairwise dart 
disjoint (due to Corollary [2]) . Applying the following lemma with k = logn and W = r then gives 
Theorem 

Lemma 7. Consider a set of objects, where each object o is assigned a positive integer weight 
w{o). Let merge{o,o') be an operation that replaces two distinct objects o and d by a new object 
whose weight is at most w{o) + w{o'). Assume that the time to execute merge{o,o') is bounded by 
0{mm{w{o),w{o')}k) for some value k. Then repeating the merge -operation on pairs of objects in 
any order until at most one object remains takes 0{kWlogW) time where W is the total weight of 
the original objects. 

Proof. We only need to consider the hard case where in beginning, all objects have weight 1 and 
at termination, exactly one object of weight W remains. Furthermore, we may assume that the 
weight of merge{o, o') is exactly w{o) + w[o') for two objects o and o' . 

Consider running the algorithm backwards: starting with one object of weight VF, repeatedly 
apply an operation split that splits an object of weight at least two into two new objects of positive 
integer weights such that the sum of weights of the two equals the weight of the original object. 
Assume that split runs in time proportional to the smaller weight of the two new objects times k. 
If we can give a bound of 0{kW \ogW) for this algorithm, we also get a bound on the algorithm 
stated in the theorem. 

The running time for the new algorithm satisfies: 

T(w) < k max iT(w') + T(w — w') + cw'} 

!<«)'< Lu>/2J 

for integer w > 1 and constant c > 0. It is easy to see that the right-hand side is maximized when 
w' = [w/2\. This gives T{W) = 0{kWlogW), as desired. □ 

Theorem 10. The region subpieces of a piece of size r can be identified in O(rlog^n) time. 

5 Adding a separating cycle to the region tree 

Above, we showed how to find a compact representation of a minimum cycle cycle C separating a 
pair of faces in a region R. This cycle should be added to the basis we are constructing and in this 
section, we show how to update the region tree T accordingly. As in the previous section, let Pr 
be the region subpiece P n i? of piece P. 

When C is added to the partial basis, R is split into two regions, Ri and i?2- Equivalently, 
in T, R will be replaced by two nodes Ri and R2. The children T R will be partitioned into 
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children Ti of Ri and T2 of i?2. Define Ri to be the region defined by the children of R that are 
contained to the left of C . Likewise define R2 to be the region for the children to the right of C. 
We describe an algorithm that finds Fi and detects whether Fi is contained by int{C) or ext{C). 
Finding J-2 is symmetric. The algorithms take 0(| J-^l log'^ n + {\Pr\ + log n) time and so we 

can identify the smaller side of the partition in 0(min{|J^i|, \ J-2\} log'^ n+(|P/j| + |9PrP) log n) time. 
This will imply that the total time for all region tree updates during the course of the algorithm is 
0(n log^n) (see proof of Lemma [t]). 

Updating the region tree Given the smaller side of the partition, we use cut-and-link operations 
to update T in 0(min{|J^i|, |J"2|}logn) additional time. See Figure [T] for an illustration. Assume, 
w.l.o.g., that Ti is the smaller set. If Fi is contained by int{C) then update T by: cutting the 
edges between R and each element in Fi, linking each element in Fi to -Ri, making R the parent 
of identifying R with R2. If Fi is contained by ext{C) then update T by: cutting the edges 
between R and each element in Fi, linking each element in Fi to a new node u, making u the 
parent of R; identifying R with Ri and u with i?2- 

5.1 Partitioning the faces 

R is represented compactly: vertices in G[R] of degree 2 are removed by merging the adjacent edges 
creating super edges. Each super edge is associated with the first and last edge on the corresponding 
path. In addition to partitioning the faces, we must find the compact representation for the two 
new regions, -Ri and i?2- 




Figure 7: C (bold cycle) is given with a counterclockwise orientation. The children of R (boundary 
given by thin cycle) adjacent and to the right of C are grey. The edges to the left of C (and not 
on C) will never reach a boundary edge of R: therefore the left of C forms int{C). Vertices of L 
are given by dark circles. 

The algorithm for finding Fi starts with an empty set and consists of three steps: 

Left root vertices Identify the set L of vertices u on C having an edge emanating to the left of 
C; also identify, for each v £ L, the two edges on C incident to v (in G, not the compact 
representation) . 

Search Start a search in R from each vertex of L avoiding edges on C or emanating to the right 
of C; for each super edge e of R visited, find the first (or last) edge e on the path represented 
by e. If there is no vertex in L take any edge e on C. 
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Add For each pair of faces /i, /2 adjacent to e, find the two children of i? in T having fi and /2 
as descendants, respectively, and add these nodes to Ti. 

This algorithm correctly builds J-i : The algorithm visits all super edges e that are strictly inside 
R and on the left side of C. Let Ai and A2 be the children of R that are added corresponding to 
edge e. Ai is a region or a face of G: let Cj be the bounding cycle. Since fi is a descendent of Ai, 
fi is contained by int{Ci). Since e is in Ci and C2: so must e. ^1 and are therefore the child 
regions of R on either side of e. 

The algorithm can easily determine if J-'i is contained by int{C) or ext[C): let /i and /2 be 
the adjacent faces of a searched edge e. Lemma [6] tells us if e is in the bounding cycle of R. If it 
is, then J-'i must be contained by ext{C): otherwise, the search could never reach an edge on the 
cycle bounding R (since edges of C are not visited) and Ti must be contained by int{C). 



Analysis 

The above-described algorithm can be implemented in 0(|J^i| log^ n + {\Pr\ + \dPji\'^)logn) time. 
Finding the left-root vertices is the trickiest part; while \L\ = 0(|J-"i|), \L\ could be much smaller 



than the number of vertices in C, even in the compact representation. We give details in Section 5.2 
Assuming that left-root vertices can be found quickly, we analyze the remaining steps. 



Search step We first need to show that the number of super edges in G[R] that have to be 
considered is bounded by \J-'i\- 

Lemma 8. The number of super edges in G[R] that has to be searched inside C is 0(|J-"i|). 

Proof. Since G has degree three, so has the compact representation of R since all faces and isometric 
cycles in G are simple. We remove all degree 2 vertices from (^[-R], so it is actually 3-regular. Hence 
the number of edges e is equal to where v is the number of vertices. By Euler's formula 

V — e + \Fi\ = 2 we get: 

le-e + \Fi\ = 2 =^ e = 3(|Fi|-2). 

□ 

The search can be done by DFS or BFS in linear time, starting with vertices of L. Given a 
super edge e found by this search, we find the first or last edge e (of G) on the path the super edge 
represents; this takes 0(1) time since e is associated with e. 

In order for this to work, we need to identify, for each v £ L, the set of super edges incident to 

V which are not on C and which do not emanate to the right of C As noted in the above proof the 
compact representation of R has degree three. It follows that no super edge incident to v emanates 
to the right of C. To identify the super edge incident to v which is not on C, we can apply Lemma[6] 
to the first edge on each of the three super edge paths. 



Add step The faces /i and /2 incident to an edge e can be found in constant time using the 
graph representation of G. Applying Lemma[6]to e (to check if J^i is contained by int{C) or ext{C)) 
takes O(logn) time. Finding the children of R that are ancestors of /i and /2 also takes O(logn) 
time using the operations jump{R, fi, 1) and jump{R, /2, 1) in the top tree for T. The total time 
spent adding is 0(|^i| logn). 
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5.2 Finding left-root vertices 

We show how to find the set L of left-root vertices along C in 0(|J^i| log'^ n + |C| logn) time where 
|C| is the number of super edges in the compact representation of C. As a result of Section [sj C 
has 0(|Pr| + |c?PrP) super edges and of three different types corresponding to: edges in exiDDG 
and iniDDG(PR), and edges and cycle paths in Pr. We will show how to use binary search to 
prune certain super edges of C that do not contain vertices of L. Assume for now that each super 
edge is on the boundary of a child region (as opposed to a child face) of R that is to the left of C. 
We will show how to relax this assumption in Section [5.2.3 

The following lemma is the key to using binary search along C: 

Lemma 9. Let P he the shortest ui-to-U2 path in G that is also a subpath of C. For i = 1,2, let Ci 
be the edge on P incident to Ui and let ri be the child-region of R that is left of C and is bounded 
by Cj. Then ri = r2 if and only if no interior vertex of P belongs to L. 

Proof. The "if" part is trivial. 

By our assumption, ri and r2 are regions, not faces. Their bounding cycles must therefore be 
isometric. If ri = r2 then Lemma [T] implies that P is a subpath of the boundary of ri: no interior 
vertex of P could belong to L in this case. This shows the "only if" part. □ 

5.2.1 Shortest path covering 

In order to use Lemma [9j we cover the left-root vertices of C with two shortest paths P and Q. 
Let r be a vertex that is the endpoint of a super edge of C. Since C is isometric, there is a unique 
edge e such that C is the union of e and two shortest paths P' and Q' between r and the endpoints 
of e. Note that e could be in the interior of a super edge of C. The paths P and Q that we use to 
cover L are prefixes of P' and Q' . 

To find e, we first find e, the super edge that contains e. Since P and Q are shortest paths and 
shortest paths are unique, the weight of each path is at most half the weight of the cycle. To find e, 
simply walk along the super edges of C and stopping when more than half the weight is traversed: 
e is the last super edge on this walk. 

Given e, we continue this walk according to the type of super edge that e is. If e corresponds to 
a cycle path, then, by definition, all the interior vertices of e have degree two in R and so cannot 
contain a left-root vertex; there is no need to continue the walk. P and Q are simply the paths 
along C from r to e's endpoints. This takes 0(|C|) time. 

If e is an edge of intDDG^P^) or extDDG, we continue the walk. We describe the process for 
miDDG(PR) as extDT>G is similar: we continue the walk started above through the subdivision 
tree of Pr that is used to find iniDDG(PR). e is given by a path of edges in the internal dense 
distance graph of Pr's children in the subdivision tree. We may assume that we have a top tree 
representation of the shortest path tree containing this path and so we can find the child super edge 
Cc that contains e using binary search taking O(log^n) time. Recursing through the subdivision 
tree finds a cycle path or edge that contains e for a total of 0(log^ n) time. 

When we are done, P and Q are paths of super edges from extDT>G or intDDG{Pji). P and Q 
each have 0(|C| -|- logn) super edges and they are found in 0(|C| + log^ n) time. 

5.2.2 Building L 

Using Lemma[9j we will decompose P into maximal subpaths Pi, . . . , P^ such that no interior vertex 
of a subpath belongs to L. Each subpath Pj will be associated with the child region of R to the 
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left of C that is bounded (partly) by Pi. We repeat this process for Q and find L in 0{k) time by 
testing the endpoints of the subpaths. 

Let e be one of the 0(|C|+logn) super edges of P. e is either an edge of exiDDG or miDDG(PR). 
Suppose e is in iniDDG(PR). We can apply Lemma [9] to the first and last edges on the path that e 
represents, and stop if there are no vertices of L in the interior of the path. Otherwise, with the top 
tree representation of the shortest path tree containing the shortest path representing e, we find 
the midpoint of this path and recurse. If e is in exiDDG, the process is similar. Adjacent subpaths 
may still need to be merged after the above process, but this can be done in time proportional to 
their number. 

How long does it take to build L? Let e be a super edge representing subpath of P and let 
m be the number of interior vertices of Pe belonging to L. Then there are m leaves in the recursion 
tree for the search applied to e. We claim that the height of the recursion tree is O(log^n). Let 
5 be some root-to-leaf path in the recursion tree. If e is in miDDG(PR), S is split into O(logn) 
subpaths, one for each level of the subdivision tree; in each level, the corresponding subpath is 
halved O(logn) times before reaching a single edge. If e is in exiDDG, the search may go root- 
wards in the subdivision tree but once we traverse down, we are in m^DDG and will thus not go 
up again. The depth of the recursion tree is still 0(log^ n). 

At each node in the recursion tree, we apply two top tree operations to check the condition in 
Lemma [9] and one top tree operation to find the midpoint of a path for a total of O(logn) time. 
The total time spent finding the m vertices of L in Q is O(mlog^n) time. If m = 0, we still need 
O(logn) time to check the condition in Lemma [9j Summing over all super edges of P, the time 
required to identify L is 0(|C| logn + |L| log^n) = 0(|C| logn + log^n), as desired. 



5.2.3 Handling faces 

We have assumed that every child of R incident to and left of C is a region, not a face. Lemma [9] is 
only true for this case: boundaries of faces need not be isometric, and so the intersection between a 
face and shortest path may have multiple components. However, notice that after the triangulation 
of the primal followed by the triangulation of the dual, every face / of G is bounded by a simple 
cycle of the form eiPie2P2GzPz where ei, 62, and 63 are edges and Pi, P2, and P3 are tiny-weight 



shortest paths (see Section 2.1). Call the six endpoints of edges 61,62,63 the corners of f. We 



associate each edge of / with the path containing it among the six paths ei, 62, 63, Pi, P2, and P3. 



We present a stronger version of Lemma 10 which implies the correctness of the left-root vertex 



finding algorithm even when children of R are faces, not regions: 

Lemma 10. Let P be the shortest ui-to-U2 path in G that is also a subpath of C. For i = 1,2, let 
6j be the edge on P incident to Ui and let ri be the child of R to the left of C and containing Ci . 
Then 

1. if neither ri nor r2 are faces, then ri = r2 if and only if no interior vertex of P belongs to L, 

2. if exactly one of ri,r2 is a face then some interior vertex of P belongs to L, 

3. if both ri and r2 are faces and ri 7^ r2 then some interior vertex of P belongs to L, 

4. if both ri and r2 are faces, ri = r2, and 61 and 62 are associated with different subpaths of ri 
then some interior vertex of P is a corner of ri or belongs to L, 

5. if both ri and r2 are faces, ri = r2, and 61 and 62 are associated with the same subpath of ri 
then no interior vertex of P belongs to L. 
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Proof. Part 1 is Lemma |9] and parts 2 and 3 are trivial. For part 4, we may assume that P is fully 
contained in the boundary of ri since otherwise, some interior vertex of P belongs to L. Since ei 
and 62 are associated with different subpaths of ri, it follows that some interior vertex of P is a 
corner of ri. For part 5, we may assume that ei 7^ 62- Then ei and 62 are on the same (tiny weight) 
shortest path in ri so P must be contained in the boundary of ri . It follows that no interior vertex 
of P belongs to L. □ 



Using Lemma 10 instead of Lemma [9j our L- finding algorithm will also identify corners of faces 
incident to P. Since each face has only 6 corners but contributes at least two vertices to L, this 
will not increase the asymptotic running time. 



5.3 Obtaining new regions 

While we have found the required partition of the children of R and updated the region tree 
accordingly, it remains to find compact representations of the new regions Ri and R2. Recall that 
we only explicitly find one side of the partition, w.l.o.g., J-i. 

To find start an initially empty graph. In the search step, we explicitly find all the super 
edges of Ri that are not on the boundary of C. Remove these edges from R and add them to 
The remaining super edges are simply subpaths of C between consecutive vertices of L. These 
edges can be added to Ri without removing them from R. 

The super edges left in R are exactly those in R2. However, there may be remaining degree- two 
vertices that should be removed by merging adjacent super edges. All such vertices, by construction, 
must be in L, and so can be removed quickly. 

That super edges are associated with the first and last edges on their respective paths is easy 
to maintain given the above construction. The entire time required to build the new compact 
representation is 0(|J^i|). 



6 Reporting min cuts 

By Theorem |6j we can report the weight of any min st-cut in constant time. We extend this to 
report the cuts themselves in time proportional to their size. Size refers to the number of edges 
in the cut. We will show how to report a minimum separating cycle for a given pair of faces in 
G* in time proportional to the number of edges in the cycle. By duality of the min cuts and min 
separating cycles, this will prove Theorem [7] 

In this section, we do not assume that the graph has degree 3. The edges added to achieve that 
may increase the number of edges in a cycle. However, we can still compute T, the region tree 
of G*, with the degree-3 assumption. We will rely only on the relationship between faces in G* , 
which did not change in the construction for the degree-3 assumption. Since cycles in the min cycle 
basis are boundaries of regions represented by T, this tree also reflects the ancestor /descendant 
relationships between cycles in the min cycle basis. 

Consider a cycle C. By Lemma [6] the set of darts in C that are not in an ancestor of C form a 
path (possibly equal to C). It follows that the set of darts in C that are not in any strict ancestor 
of G also form a path, denoted P{G). Using the next lemma, we can succinctly represent any cycle 
using these paths. See Figure [8] for an illustration. 

Lemma 11. Let G = Gq, Gi, G2, . . . be the ancestral path to the root ofT for G . G can he written 
as the concatenation of the path P{G), prefixes of P{Gi),P{G2), ■ ■ ■ ,P{Gk-i), a subpath of P[Gk) 
and suffixes of P{Gk-i), ■ ■ ■ , P{G2), P{Gi), in that order. 



20 




d^(C^) d2(C^) 



Figure 8: A succinct representation of a cycle in the min cycle basis. 



Proof. Let be the root-most cycle that shares a dart d with C: d is in P(Cfc). By Lemma |6| the 
intersection of C with Ck is a single path: it must be a subpath P' of P{Ck)- Let Q = C \ P' . By 
the definition of P{Ck~i), the start of Q must be the start of P{Ck~i) and the end of Q must be 
the end of P{Ck-i)- The remainder of the proof follows with a simple induction. □ 

Let di{C) be the last dart of C before P{C) and let d2{C) be the first dart of C after P{C). 
Suppose additionally that we know, for every dart d, the cycle C{d) for which d G P{C). 



6.1 Finding the min separating cycle 

If we are given di{C), d2{C), and P{C) for every cycle (node) in T and C{d) for every dart d, we 
can find the minimum /^-separating cycle C in 0(|C|) time by the following procedure. First we 
can find the node in T corresponding to C in 0(1) time using the oracle (Theorem [6]). To find C, 
walk along C, starting with P{C), until you reach the end. Let Ci = C{d2{C)) and walk along 
P{Ci) starting with d2{C). Suppose we are at dart d along C. Let Ci = C{d). Walk along P{C.i) 
until either you reach its end or you hit dart di[Ci^i). In the first case, continue the process with 



d2{Ci). In the second case, continue the process with the first dart of P(Cj_i). By Lemma 11, this 
process will eventually reach the start of P{C). 



6.2 Preprocessing step 

It remains to show how to precompute (ii(C), d2{C), P{C) for every cycle in T and C{d) for every 
dart. We find these using the top tree representation of T with 0(n log n) preprocessing time. 

Let fn and fr be the faces to the left and right of a dart d. Then it follows easily from Lemma [6] 
and the clockwise orientation C{d) that C{d) is the bounding cycle of region ju'mp{lca{fi, fr), fr, 1) 
and can be found in O(logn) time. 

We can easily construct P{C) from the set of darts with C{d) = C. The ordering can be found 
just using the endpoints of these darts so that we can walk along P{C) as required in the previous 
section. 

To find di(C) and 1^2 (C), we work from leaf to root in T as in the bottom-up algorithm of 



Section 4.3, We will show how to find d2{C). Finding di{C) is symmetric. For cycle C we 
can easily find the last dart di of P{C). Consider the darts do leaving the endpoint of P{C) in 
counterclockwise order, in the embedding, starting with the reverse of dg, we test if do is on the 
boundary of C using Lemma |6] in O(logn) time. As we test darts we remove them from further 
consideration as they will be in the interior of all ancestor cycles. In total, this takes 0(n log n) 
time. 
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7 Handling Holes and Disconnected Pieces 



In Section [l.2| we introduced the notion of recursive subdivision of the planar graph. We made two 
simplifying assumptions, namely that each piece is connected and does not contain holes. In the 
following, we show how to avoid these assumptions. 

7.1 Disconnected Pieces 

The problem with a disconnected piece P is that a minimum separating cycle for a pair of faces 
touching two distinct connected components of P need not share any vertex with P. In such a 
case, the algorithm in Section [3] cannot be applied to find this cycle. 

To deal with this, we first modify the definition of the recursive subdivision by allowing a piece 
to have more than two children. More precisely, if a piece is split into k connected components 
by the application of Miller's cycle separator theorem, we assign k children to it in the recursive 
subdivision, one for each connected component. This ensures that pieces are connected. Also, the 
size of pieces will still go down geometrically when moving down the recursive subdivision tree so 
the height of the tree remains O(logn). Furthermore, intDDG can be found in O(nlog^n) time 
using the algorithm of Fakcharoenphol and Rao. 

In order to compute exiDDG within the same time bound, we need to generalize the top-down 
algorithm contained in (the proof of) Theorem [2] Consider a piece P and suppose extDDG{P) 
has been computed. Let Pi,. . . ,Pk be the children of P in the recursive subdivision. For each 
i, we compute extDT)G{Pi) with the Dijkstra variant of Fakcharoenphol and Rao on the union of 
extDDG{P) and intDDG{Pj), j = l,...,k, j / i. We apply this algorithm for each boundary 
vertex of Pi as source, giving all the pairwise distances needed in extDT)G{Pi). Since the total 
number of boundary vertices in Pi,...,Pk is 0{y^\P\), each application of the algorithm runs 
in 0{y^\P\ log^ n) time. Summing over all boundary vertices of all children of P, the total time 
to compute extDBGiPi), extDDG{Pk) is 0{\P\ log^ n). Hence, extDBG can be computed in 
0(nlog^ n) time. 

Now, for the algorithm in Section [3j we proceed as before, separating pairs of faces for pieces in 
a bottom-up fashion in the recursive subdivision tree. Consider a piece P with k children Pi, . . . ,P^. 
and assume that the algorithm has been applied to each of them. We let 6 Pi , ■ ■ ■ , SPk denote the 
boundary vertices of Pi, . . . , P^, respectively. 

To separate unseparated faces touching P, let us first consider the simpler problem of separating 
faces touching two distinct children, say. Pi and P2. We find region subpieces for their union with 
the algorithm in Section |4| This takes 0((|Pi| -|- |P2|)log^n) time. A proof similar to that of 
Lemma |4] shows that the number of region subpieces is Od^Pil -|- |(5P2|)- 

For each such region subpiece Pr, there is exactly one pair of faces that needs to be separated, 
one touching Pi and one touching P2. Since P is connected, we know that a minimum separating 
cycle for this pair must intersect P. Hence, this cycle can be represented as the union of edges from 
Pr, edges from exiDDG(P), and/or edges from miDDG(P3), . . . , intDDG(Pfc). 

To find this cycle, we first find (a compact representation of) the shortest path X between 
the two faces as in Section [sj This takes 0(|PR|logn -|- \/\P\ log^ n) time using a combination 
of normal Dijkstra inside Pr and the Dijkstra variant of Fakcharoenphol and Rao in the dense 
distance graphs. After splitting Pr and the dense distance graphs along X as in Section |3.2[ we 
obtain cut-open graph Pji^ from Pr, intDDGx{Pj) from intDDG{Pj) for each piece Pj, j / 1,2, 
and extDDGx(P) from exiDDG(P). In the following, let H be the union of all these resulting 
graphs. 

We use a variant of Reif 's algorithm to find a minimum separating cycle. Let x be the midpoint 
of (the compact representation of) X. We find a minimum separating cycle in H that crosses X 
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at X, using a combination of normal Dijkstra inside Pr and the Dijkstra variant of Fakcharoenphol 
and Rao outside Pr. Then we recurse on the two subgraphs. This takes a total of 0(|Pr| log^ n + 
^J\P\\og^ n) time. 

Each subgraph generated by this algorithm now contains exactly one edge of X. Had all edges 
of X been edges of G, we would have found a minimum separating cycle. However, some edges of 
X may belong to dense distance graphs and may thus represent paths in G so additional separating 
cycles need to be found. In the following, we consider a subgraph of H containing an edge 
e G X from intDT)Gx{Pj), j 7^ 1,2. The case where the edge belongs to extDDGx{P) is handled 
in a similar way. 

We wish to find a separating cycle that crosses X somewhere on (the path in G represented 
by) e. We cannot use the idea from Section 3.4 of looking at all pairs of boundary vertices of 
intDT)Gx{Pj) since the number of pairs that need to be considered may be much larger than the 
number of pairs in Pr. 




Figure 9: Finding an /(7-separating cycle through a super edge e of X that belongs to a sibling 
piece Pj (j / 1,2). e splits (a subset of) the boundary vertices of Pj into two sets: Ui is given by 
the solid vertices and C/2 by the hollow vertices. A separating cycle through the middle vertex m 
of Ui (according the order along dUi) is given by the dotted cycle. 



Instead, we will apply divide- and-conquer on boundary vertices in intDDGxiPj)- Consider 
the set Ui x U2 of all pairs {ui,U2) of boundary vertices in intDDGx{Pj) such that edge (ui,U2) 
crosses X exactly once and does so in e, see Figure [9| In the ordering of boundary vertices induced 
by the external face of Pj, we pick ui as the middle vertex in Ui. We then compute a shortest 
non-empty path Cui from ui to itself in subgraph such that the first edge on the path is of the 
form ei = (ui,«2), «2 G f^2- This path defines a separating cycle since it crosses X exactly once. 

For divide-and-conquer to work, we need to separate our problem in two. To see that it 
does, suppose a minimum separating cycle in crosses Cm- This cycle can also be decomposed 
into an edge e'l = {u'i,U2) G Ui x U2 and a shortest path in He from u'2 to u'l, we denote it by 
C^. It must cross in ei since otherwise it would have to cross G^ \ {ei} twice, contradicting 
the shortest path uniqueness assumption. In fact, the same argument shows that C„/ must cross 



Gui exactly once in ei and exactly once in Cui \ {ei}, see Figure 10 But then either we can find 



a separating cycle through ui which is shorter than Cu^ or a shorter separating cycle than C^/ . It 
follows that Gui separates our problem, as desired. 

We also need to show how to efficiently find Gui- Take a copy u'l of ui and for each U2 G U2, 
add edge {u'i,U2) with weight equal to that of (ui,U2) in intDDGxiPj)- Then we obtain as a 
shortest path from u'l to Ui. Running time is 0(|Pr| log^ n + \/\P\ log^ n). 

Having found G^ , we recurse on the two subgraphs of He ■ Clearly, a minimum separating cycle 
crossing e can be realized as one such cycle Gui- Since we always pick the midpoint of Ui, the 
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Figure 10: Cycle C„'^ (grey) cannot cross Cui (dashed) by uniqueness of shortest paths: the dashed 
grey cycle would be shorter. Recall that e, ei and e'^ are all super edges. 

recursion depth is logarithmic in |J7i|. It follows from this and the above that the total time to find 
a minimum separating cycle crossing X exactly once is 0{\Pr \ log^ n + -^/|P| log^ n). 

We apply the above algorithm for all region subpicccs Pr. This takes 0((|Pi| + |P2|)logn + 
{\SPi\ + \6P2\)y/\P\log^ n) time. This process is repeated for pairs (P3, P4), (P5, Pq)! • • • using a 
total of 0(|P|logn+ y|P[y|P| log^ ra) = 0(|P|log^n) time. 

Now, each graph Pi U P2,P3 U P4, . . . contains no unseparated pairs of faces. We can then 
perform another run of the above algorithm but for pieces P[ := Pi U P2, P2 := P3 U P4, . . . instead 
of Pi, P2, . . . , Pfc. For each run, we halve the number of pieces so after O(logn) runs, we have 
separated all faces touching P. Each run takes 0(|P| log^ n) time for a total of 0(|P| log"^ n) time. 
Over all pieces P in the recursive subdivision, this is O(nlog^n). 

7.2 Holes 

Above, we showed how to deal with disconnected pieces. We now show how to adapt our algorithm 
to handle holes. As shown by Fakcharoenphol and Rao, we may assume that the number of holes 
in a piece is bounded by a constant. 

Computing intDDG can still be done in 0{n log^ n) time since the algorithm of Fakcharoenphol 
and Rao also works in the presence of holes. To compute extDDG within the same time bound, we 
will modify the above algorithm dealing with disconnected pieces. Consider a (connected) piece P 
with h holes and suppose exiDDG(P) has been computed. Note that exiDDG(P) is the union of 
h + 1 graphs: one graph defines pairwise distances for boundary vertices on the external face of P 
and each of the other h graphs defines pairwise distances for boundary vertices in a hole. 

Let S be the separator cycle applied to P to give its (connected) children Pi, . . . ,Pfc in the 
construction of the recursive subdivision of G. We can divide these children into two types, those 
contained in the interior and those contained in the exterior of S. Let Vint resp. Vext be the set of 
children of the former resp. latter type. 

Consider a child Pj € Vint- If Pi contains a hole, it must be a hole of P. Hence, we only need to 
worry about the external face of Pj when computing exiDDG(Pj). All pairwise distances between 
boundary vertices on this face can be obtained with the Dijkstra variant of Fakcharoenphol and 
Rao on the union of the internal dense distance graphs of the siblings of Pj and the subgraph of 
extDDG{P) that excludes the holes of Pj. 

Now, consider a child Pj G Vext- We compute pairwise distances for boundary vertices on 
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the external face of Pi as above. A hole of Pi is either a hole of P or the new hole H (if any) 
introduced when separating P with separator cycle S. We only need to worry about H. The 
pairwise distances between boundary vertices on H can be obtained from the union of the internal 
dense distance graphs of siblings of Pi belonging to Vint and the subgraph of exiDDG(P) defined 
by the holes of P contained in the interior of H. 

It is easy to see that the above algorithm computes extDDG in O(nlog'^n) time. Now, let us 
consider our algorithm to find region subpieces and to separate faces in Sections [3] and [4j It is 
not hard to modify it for holes. Each hole can be handled in a way similar to how we handle the 
external face (simply consider an embedding of the graph in which the hole is the external face 
whereas the original external face becomes a hole). The only part that needs focus is the fact that 
the shortest path X, used in Section [Sj can go through several holes. For each of these holes we 
need to find external cycles in Section [3 . 4| separately. However, the total number of boundary nodes 
in all holes plus the external face is equal to the total number of boundary nodes in a piece. Hence, 
the total running time of all searches does not increase. It follows that our bounds also hold in the 
presence of holes. 



8 Lex-shortest paths 

Let 1(7 :£'—)■ M be the weight function on the edges of G. So far, we have assumed uniqueness of 
shortest paths in G between any two vertices w.r.t. w. We now show how to avoid this assumption. 
We assume in the following that the vertices of G are given indices from 1 to n. For V' C V, define 
I{V') as the smallest index of vertices in V' . 

As shown in [TJ, there is another weight function w' on the edges of G such that for any pair 
of vertices in G, there is a unique shortest path between them w.r.t. w' and this path is also a 
shortest path w.r.t. w. Furthermore, for two paths P and P' between the same pair of vertices in 
G, w'{P) < w'{P) exactly when one of the following three conditions is satisfied: 

1. P is strictly shorter than P' w.r.t. w, 

2. P and P' have the same weight w.r.t. w and P contains fewer edges than P' , 

3. P and P' have the same weight w.r.t. w and the same number of edges and I(y{P)\V{P')) < 
I{V{P')\V{P)). 

A shortest path w.r.t. w' is called a lex-shortest path and a shortest path tree w.r.t. w' is called a 
lex-shortest path tree. 

The properties of w' allow us to apply this function instead of w in our algorithm. In the 
following, we show how to do this efficiently. 

We first use a small trick from [7]: for function w, a sufficiently small e > is added to the 
weight of every edge. This allows us to disregard the second condition above. When comparing 
weights of paths, we may treat e symbolically so we do not need to worry about precision issues. 
The tricky part is efficiently testing the third condition. 

We need to make modifications to every part of our algorithm in which the weights of two 
shortest paths are compared. All such comparisons occur when we (1) apply the Dijkstra variant 
of Fakcharoenphol and Rao ^ and (2) find a shortest path covering of an isometric cycle C in 
Section 15.2.11 



Note that in making the graph degree three (Section 2.1 ), we used an e weight on the introduced 
edges. To combine this with lex-shortest paths, we make e 3> e so that shortest paths will still 
be chosen to have the fewest edges (in the original graph) and run along the trees introduced to 
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make vertices degree three. Using both e and e in the symboHc comparisons does not increase the 
asymptotic run time. 



8.1 Dijkstra 

Let us first adapt Dijkstra to compute lex-shortest paths. The type of shortest path weight com- 
parisons in [5] are of the form D{u) + d{u, v) < D{u') + d{u' , v), where u, v, u' , and v' are vertices, 
D{u) and D{u') are the distances from the root of the partiahy built tree to u and u' , respectively, 
and d{u,v) and d{u',v) are the weights of edges {u,v) and {u',v). Here, an edge can belong to G, 



be a cycle edge (see Section 4.3), or an edge of an external or internal dense distance graph. 

For simplicity, assume first that all edges considered by Dijkstra belong to G. Then we can test 
whether D{u) + d{u, v) < D{u') + d{u' ,v) as follows. Let T be the partially built shortest path tree 
rooted at a vertex r and let Q and Q' be the lex-shortest paths from r to n and n', respectively. 
If the first two lex-shortest conditions are inconclusive, we need to check if I{V{Q) \ V{Q')) < 
IiV{Q')\V{Q)). 

Letting a = lca{u,u',) in T, V{Q)\V{Q') is the set of vertices on Q[a,n], excluding a. Similarly, 
V{Q')\V{Q) is the set of vertices on Q'[a, u'], excluding a. It follows from this that by representing 
T as a top tree, we can find the smallest index in the two sets in logarithmic time. Using top trees, 
we can also deal with cycle edges in the same way by keeping, for each such edge, the smallest index 
of its interior vertices. These indices can be found during the construction of region subpieces in 



Section 4.1 without an increase in running time. 



8.2 Internal dense distances 

We also need to handle edges from internal and external dense distance graphs. Let us first 
consider the problem of computing lex-shortest path trees in iniDDG. As before, we compute 
shortest path trees for pieces bottom-up. Let P be a piece with children Pi and P2 and assume 
that we have computed lex-shortest path trees in both of them. Assume also that every edge in 
miDDG(Pi) U intDT)G{P2) is associated with the smallest index of interior vertices on the path in 
G that the edge represents. This information can be computed bottom-up during the construction 
of intDDG without increasing running time. 

Let T be a partially built shortest path tree in P and with the above defintions, consider the 
problem of testing whether D{u)+d{u, v) < D{u')+d{u' , v). Let (a, Ua) and (a, u'^) be the first edges 
on Q[a, u] and Q'[a, u'], respectively. Define Qg and Q'q as the paths in G represented by Q and Q' , 
respectively. Let h = I{V{QG[ua,u])), i[ = I{V{Q'cK,u'])), i2 = I{V{QG[a,Ua]) \ V{Q'c[a,u',])), 
and i2 = I{V{Q'Q[a,u'^]) \ V{QG[a,Ua]))- By definition of a, QG[ua,u] and Q'q[u'^,u'] are vertex- 
disjoint. Hence, we need to find the four indices and check if minjii, 12} < min{z'^, ig}. 

Each edge of T belongs to miDDG(Pi) U intT)DG{P2) and is thus associated with the smallest 
index of interior vertices on the path in G represented by that the edge. Top tree operations on T 
as above then allow us to find ii and i'^ in logarithmic time. 

To find 12 and we consider two cases: (o, Ua) and (a, u'^) belong to the internal dense distance 
graph for the same child of P or they belong to different graphs. In the first case, assume that, 
say, (a, Ua), {a,u'g^) G miDDG(Pi). Then we can decompose these two edges into shortest paths in 
the same shortest path tree in miDDG(Pi) and we can recursively find ^2 and In the second 
case, assume that, say, {a,Ua) S zniDDG(Pi) and {a,u'^) S intT)DG{P2). Since the lex-shortest 
paths representing these edges in intDDG(Pi) and intDDG{P2) are edge-disjoint and since T is a 
partially built lex-shortest path tree in P, Qg[o--, Ua\ and it^] share no vertices except a. Thus, 

12 is the smallest index of vertices in V{QG[o-,Ua\) \ {a} and we can obtain this index in constant 
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time from the index of Ua and the index associated with edge (o, Ua) which is the smaUest index of 
interior vertices on Qda^Ua]. Similarly, we can find i'2 in constant time. 

Since the subdivision tree has O(logn) height, the recursion depth of the above algorithm is 
O(logn), implying that we can determine whether D{u) + d{u,v) < D{u') + d{u',v) in O(log^n) 
time. Hence, lex-shortest path trees in iniDDG can be computed in a total of O(nlog^n) time. 



8.3 External dense distances 

Computing lex-shortest path trees in ex^DDG within the same time bound is very similar so we 
only highlight the differences. Having computed lex-shortest path trees in miDDG bottom-up, we 
compute lex-shortest path trees in extDDG top-down. For a piece P, we obtain lex-shortest path 
trees from lex-shortest path trees in its sibling and parent pieces. We can then use an algorithm 
similar to the one above to find lex-shortest path trees in P. At each recursive step, we either go 
up one level in exiDDG or go to intDDG. It follows that the recursive depth is still O(logn) so 
lex-shortest path trees in extDDG can be found in 0(n log^ n) time. 



8.4 Distances in region subpieces 

We also apply Dijkstra in Section [3] for a region subpiece Pr. First, we computed shortest path X 
between two faces of the region subpiece in 0((|5PrP + |Pr|) log |-Pr|) time using Dijkstra. With 
an algorithm similar to the one above, we can instead compute a lex-shortest path between the 
two faces in 0((|5Pijp + |Pr|) log \Pr\ log^ n) time. Next, we cut open the region subpiece along 
this path and computed a shortest path in Gx between each pair of duplicate vertices on the two 
copies of X. Using Reif 's algorithm, we obtained a time bound of 0((|9Pr|2 + |P^|) log^ \Pr\). We 
need to find similar shortest paths but with respect to weight function w' . 

For each pair of duplicate vertices corresponding to a vertex x £ X, we assign the index of 
X to both of them. Then the lex-shortest path between the two vertices in Gx w.r.t. w must be 
the shortest path w.r.t. w' . Hence, using the above Dijkstra algorithm, the total time spent on 
computing shortest paths w.r.t. w' for Pr is 0((|Pr| + log^dP^I) log^ n). 



8.5 Shortest path coverings 



In Section 5.2.1 , we gave an algorithm to find the unique edge e = (u, v) on isometric cycle C such 
that the two shortest paths from a fixed vertex r on C to u and to v cover all vertices of C and 
all edges except e. We showed how to do this in 0(|C| + log^n) time, where |C| is the size of the 
compact representation of C obtained in Section |3j We need to modify the algorithm to deal with 
lex-shortest paths. 

Recall that to find e, a linear search of the super edges of C from r was first applied to find the 
super edge e of C such that the shortest path in G representing e contains e. As above, we may 
assume that every super edge of C is associated with the smallest index of interior vertices on the 
path it represents. Hence, by keeping track of the smallest interior vertex index for super edges 
visited so far in the linear search as well as the smallest interior vertex index for edges yet to be 
visited, we can find e in 0(|C|) time w.r.t. lex-shortest paths. 

Having found e, we need to apply binary search on a path representing e in a lex-shortest path 



tree. We do this by first finding the midpoint of this path as in Section 5.2.1 If the two halves 
have the same weight and the same number of edges, we can use a top tree operation on each 
half to determine which half has the smallest index. It follows that all binary searches to find e 
take O(log^n) time. The total time to find e is thus 0(|C| + log^n), which matches the time in 
Section 15.2.11 
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As a result, we have: 



Theorem 11. The Theorems^ through^ hold without the shortest-path uniqueness assumption 
with only an additional 0(log^ n) factor in the preprocessing time. 

9 Concluding remarks 

We gave an oracle for min st-cut weight queries in a planar undirected graph with non-negative 
edge weights. Construction time is O(nlog^n) and space requirement is 0(n log n), where n is the 
size of the graph. The actual cut can be reported in time proportional to its size. We obtained 
this oracle from a Gomory-Hu tree algorithm with the same time and space bounds. Previously, no 
subquadratic time algorithm was known. We also showed how to compute an implicit representation 
of a minimum cycle basis in 0(n log^n) time and O(nlogn) space and an explicit representation 
with additional 0(C) time and space where C is the size of the basis. 

In order to obtain our results, we needed shortest paths to be unique. We showed how to 
deterministically remove this assumption at the cost of an extra log^ n-factor in running time. 
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